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Abstract 



This note describes five subjects of some interest for the density func- 
tional theory in nuclear physics. These are, respectively, i) the need for 
concave functionals, ii) the nature of the Kohn-Sham potential for the 
radial density theory, iii) a proper implementation of a density functional 
for an "intrinsic" rotational density, iv) the possible existence of a poten- 
tial driving the square root of the density, and v) the existence of many 
models where a density functional can be explicitly constructed. 

> 

, Preliminary considerations : The nuclear Hamiltonian used here reads, 



i i>j \ i I 



(1) 



will be clar- 



ON . 

where 7*j and pi are the coordinates and momenta of identical nucleons, with 
spin and isospin labels understood unless necessary, m is their mass, is the 
two-body interaction, Galilean invariant. A term, J2i>j>k w ijki f° r three-body 
' interactions, can be added if needed, provided it is also Galilean invariant. The 

role of the additional terms, a center-of mass trap with strength e, and particle 

. \a B~ 

number fluctuation terms, with a positive definite matrix, 

LP 7. 

ified in time in the following. The operators N and Z are the neutron and pro- 
ton number operators. Since they commute with the physical part of H, namely 
Hphys = J^Pi v iji their presence leaves the eigenstates intact and only modi- 
fies the spectrum by trivial quantities, which can be easily subtracted eventually. 
Similarly, the center-of-mass trap commutes with that internal dynamics gov- 
erned by the internal Hamiltonian, Hi n t — ^2 Pi / (2 m ) ~ (^2 Pi) 2 /(2-Am)+2 Vij ;, 
where A = N + Z is the mass number. The internal dynamics is not perturbed, 
the center-of-mass motion factorizes out and the corresponding zero-point en- 
ergy can be trivially subtracted eventually. 
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1 Need for concavity 



At the core of density functional (DF) theory, there is the process, analyzed 
by Levy and by Lieb [1] , of a minimization of energy under the constraint of a 
given density, 



Here T> is any physical density operator in the Fock many-body space, the 
nuclear Hamiltonian H is second quantized if necessary, the trace Tr is taken in 
the many-body space, the arrows => are shorthand for the constraints, 



where ai and ap are the usual creation and annihilation operators at position 
r, with spin and isospin labels understood. It is also understood that p has two 
components, a neutron and a proton one. Spins are most often summed upon. 
The minimization, Min, may be understood as an infimum, Inf, if necessary, 
depending on fine details of the considered density operators, projectors V = 
\^){^\ or, in obvious notations, ensemble mixtures V = ^2 n w n \^f n )(^f n \. Fine 
details of the representability of p are also understood in this paper. 

It is well known, from the Hohenberg-Kohn formulation [SJ, that the con- 
straint, T> => p, can be implemented by means of auxiliary potentials u n ,u p , 
with Lagrange terms, J dru n (r) p n (r), J dfu p {r) p p (r), and subsequent Legen- 
dre transforms. 

Before the true spectrum of H exhibits simultaneous eigenstates of N, Z 
and Hi n t, nothing prevents the theory from considering density operators T> 
that are not eigenstates of the particle number operators. It is therefore also 
useful to consider, for instance, pairing densities for opposite spin particles, 
k(t) = Tr at + at _ T) , and minimize the energy under both the p and the k 
constraints. This generalizes functionals F[p] into functionals F\p, n]. See [3] 
for specificities associated with zero range pairing. 

A mass formula is concave if any second difference such as, E(N + 1, Z) — 
2 E(N, Z) + E(N -1,Z), E(N, Z + l)-2 E(N, Z) + E(N, Z-l), and similarly 
with any range in any direction of the N, Z plane, is positive definite. It is 
a fact of life, however, that the table of nuclear ground state energies is far 
from providing systematically positive second differences. Not only there is a 
significant staggering effect in binding when comparing neighboring odd and 
even nuclei, because of pairing, but even more severe deviations of concavity 
occur almost every time a magic shell or subshell closes. 

This lack of concavity has a severe consequence for a mass formula such 
as that provided by Eq. {2}. Assume indeed that, in the neighborhood of 
an N, Z nucleus, there are two integers n, z so that the experimental second 
difference of ground state energies, Em+u,z+z — 2 Sjv.Z + Ex-n.z-z is negative. 
Let ^N+n.z+z and ^N-n.z~z be corresponding ground state wave functions. 



E(N, Z) = Min p=>7V , z F[p] , F[p] = Min v=>p 



Tr H V. 



(2) 




(3) 
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(If ground states degenerate because of a spin different from zero, any magnetic 
label will do.) Then the ensemble mixture, 



provided the eigenvalues of 



are large enough. A rough inspection [3] 



V = \ {\^N+n,Z+z)(^N+n,Z+z\ + |* N-n,Z-z) (* N-n,Z-z\) , 

induces a density, \ {pN+n,z+z + PN-n,z-z) > with correct average particle num- 
bers TV and Z. But it returns an absurd energy, i (E^ +n z+ z + EN—n,z—z) > 
lower than z- 

In Eq. fT]). the terms quadratic in N and Z, however, do create concavity, 

'a f3~ 

7. 

of the table of nuclear bindings indicates that eigenvalues of a few MeV at most 
are enough to ensure concavity. 

One might argue that the Kohn-Sham (KS) determinant [5] defines particle 
numbers well, but the argument does not extend to DF calculations with pairing. 
Indeed, procedures for particle number "projection" had to be designed [6] [7|. 
It might be that energy density functional (EDF) theory 8J, with its several 
quasi local densities, is protected against such a syndrome related to ensemble 
densities, because EDFT carries more information about the nucleus, but the 
question of concavity, up to my knowledge, has not been discussed explicitly. 
Nor has it been raised in the "no Hamiltonian" context. 

It is clear that the "quadratic N, Z terms" in Eq. (fTJ) are there to penalize 
particle number fluctuations. They relocate energy minima at strict eigenval- 
ues of N and Z. This seems to make an argument to advocate Hamiltonian 
approaches to DFT and be cautious about non-Hamiltonian ones, incidentally. 
Anyhow a problem remains : how does one implement such terms in the con- 
struction of a DF? It is likely that a formalism including pairing correlations n, 
sensitive to the presence of such terms, will be needed. The precise adjustment 
of DF and EDF theories to rigorously ensure concavity is an open problem. 



2 Radial Kohn-Sham theory 

Whether one considers the usual Hamiltonian H p h ysi or the internal one, Hi nt , 
or that full H, see Eq. (jTJ) , which is adapted to DFT, such operators are scalars 
under rotation. Accordingly, the operation, TrHV, which defines the energy, 
see Eq. J2j, is sensitive to only the scalar part of T>. (Notice that the normal- 
ization, TrP = 1, also selects the scalar part of T>.) The energy minimization 
can, therefore, be restricted to a subset of purely scalar 2?s, without any loss 
of information [9]. For nuclear ground states that have a non zero spin </, the 
degenerate magnetic multiplet can be combined into a scalar density operator, 
the mixture T>o — J2m \ JM){JM\/(2J + 1), with, obviously, the ground state 
energy given by, Ej = Tr Hi nt T>q . The same mixture results from DFT at finite 
temperature [10j . at the zero temperature limit. 

Since the energy minimization occurs within a subspace of scalar density 
operators, the corresponding one-body densities p(f) are just radial profiles, 
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Po{r). The final result is a scalar functional, F n [p n }. Hence an open problem : 
for such monopolar profiles po, what is the signature for deformation? It is likely 
that the slope of po at the nuclear surface will be weaker for deformed nuclei 
than for spherical ones, but precise criteria distinguishing hard from soft nuclei 
are desirable, not to mention criteria identifying triaxial nuclei, halo nuclei, etc. 
Can one use properties of moments of the profiles? Such a theory with radial 
profiles, in one dimension rather than the three dimensions of usual nuclear 
(E)DF theories, provides a considerable simplification of numerics, but demands 
somewhat subtle criteria for the identification of deformations. 

Another open problem consists in adapting the KS procedure to this radial 
world. Slater determinants, made of LS orbitals, \nlma), or jj ones, \nljm), 
orbitals provided by scalar mean fields uo(r), are never scalar wave functions, 
except at subshell closure. It seems therefore necessary to redefine the KS kinetic 
DF from a formula similar to Eq. ((2]), 



Fkso[po] = Minx>' = >p Tr 



E^/(2m) 



K (4) 



where V Q means scalar mixtures of determinants made of the LS or jj orbitals. 
For instance, given a neutron number N, assume a filled neutron core accom- 
modating N' neutrons and an open subshell with spin j to accommodate the 
(N — N') neutrons left. Then T>' can represent a mixture, with equal weights, 
of all combinations of (N — N') orbitals within the (2j + 1) available ones. For 
a justification of the equal filling approximation, see [IT] . 
Then a KS potential ukso can be defined as, 

( v S(F - F KSQ ) 

ru K so(r) = j— 7-7 , 5 

rdp Q {r) 

and can be further simplified if a Hartree subtraction makes sense for nuclear 
forces. Finally a spherical set of orbitals must be generated, from the potential 
v-ksq, to obtain self-consistency of po- 

Little is known at present in nuclear physics about this redefined kinetic DF, 
Eq. (|4"j). It is time to import, and maybe readjust, the results that are known 
in atomic physics [T2] [T3]. While the early literature of nuclear physics is not 
devoid of effective spherical potentials for shell models, and while there has been 
recently interesting attempts to describe spherical nuclei with empirical EDFs, 
see for instance [14] , very little is known about ukso either. It may be also worth 
comparing such "spherical" EDFs with those used for deformed nuclei, see for 
instance |15j . The program of a radial KS theory is a largely open problem. 



3 Functionals for intrinsic rotational states 

Let us first avoid a frequent semantic confusion between "intrinsic" and "inter- 
nal" . In this note, the word "intrinsic" shall be used for a state which has to 
be simple, and out of which rotational states can be obtained. Simplicity means, 
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in particular, that one uses microscopic degrees of freedom, in the laboratory 
frame. The state $ is typically a Slater determinant or a Hartree-Bogoliubov 
state. Rotational states are calculated from the usual algebra with rotation op- 
erators and rotation matrices. Skipping all technical details, the main formula 
in this algebra reads, in transparent notations, 



Hik= dtl>d6dpDi IK (^,6,y>) ex.p(-ii;J z )exp(-i6J y )exp(-i<pJ z )\<f>). (6) 



(If necessary, a subsequent configuration mixing, * = J2k 
prove the calculation of a physical state.) The main point of interest is that 
the states &mki because of their quantum numbers, contain correlations that 
$ must not contain in order to be "simple" . 

It is well known that, except at the limit of very hard rotators, there is 
no decoupling possible between the three rotation degrees of freedom and any 
residual set of (3A — 3) "internal" degrees of freedom describing any "internal" 
motion inside the rotator. Such is not the case for the translation group, where 
the internal Hamiltonian, H int = J2Pi/(^ m ) ~ (J^Pi) 2 /{2Am) + J2 v ij> an d 
the trapped center-of-mass one, H com = {J^Pi) 2 /(2-Am) + e (J2^i) 2 > drive 
separate dynamics. The subtle relation between internal density and density in 
the laboratory and the theories for the corresponding DFs make a closed problem 
[16], which is, therefore, out of the scope of this note. The word "internal" shall 
be kept for the context of the translation group and the word "intrinsic" shall 
be used here, specifically, in the context of the rotational group. 

Assume that every & MK can be a good approximation to an eigenstate of 
H and consider the two-step minimization, 



where p is the density of 4>, and <& is restricted to a class of simple states. This 
defines a density functional theory [17j in terms of the intrinsic density, not the 
density of that (hopefully physical) state $ J MK which provides the energy. But, 
obviously, this DF, Fjmk[p], depends on that subset to which $ is restricted. 
Moreover, it depends on J, M , K. 

There are not many subsets available for Hartree-Bogoliubov solutions do 
not make a bad first choice. But then the essentially desirable virtue of <E> would 
be to be practically independent of J, M, K. Can this be true? Hence the open 
problem : under which conditions can Fjmk define just one density p valid for 
a large number of combinations JMK7 How does this intrinsic DF relate [18] 
to the traditional three-dimensional (E)DFs used for deformed nuclei? 

4 About the square root of the density 

A theorem, studied by many authors in atomic or molecular physics, see in 
particular [19], states that, given the density p(r) of an eigenstate of H with 



Ejmk = Min p Fjmk [p] , F JMK [p] = Min$ =>p 



«k\H\<k) 
(KkHik) 



(7) 
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energy Ea, there exists a local potential, u sqr , that drives the square root of p, 



l-h 2 A P /(2m) + u sqr {r)\ ^fpj¥) = (E A - E A -x) ^/W) ■ 



(8) 



Here Ea-i is the ground state energy of the system obtained by removing one 
particle. (When Ea corresponds to the ground state, the difference, (Ea — 
Ea-i), is a separation energy.) The details of the proof, found in [T5], will not 
be discussed here. The theorem was recently generalized for mixtures of degen- 
erate states [20] . Let us only mention that part of u sqr consists in a Hartree-like 
convolution of the two-body interaction with a density in the (A — 1) space. 
The theorem was established for electrons in atoms or molecules, where Vij is 
the Coulomb repulsion, local. Locality plays a role in the proof of the theorem. 
It seems, therefore, that an extension of the theorem with nuclear interactions, 
which can be seriously non local, is excluded. But the nuclear DFT makes an 
extensive use of locality or quasi-locality. The KS potential, in particular, is 
a functional derivative, local. Moreover, effective forces can be local, optical 
potentials can be local, Hartree potentials can be local, up to reasonable ap- 
proximations at least. Can there be a generic form of an approximate, local 
u sqr , valid for many nuclei? Would it relate to the nuclear DFT? 

5 Constructive models and their polynomial zoo 

Everyone knows that existence theorems for DFs do not provide simple, explicit 
constructions [21] [22j . However, a direct, constructive approach has recently 
been proposed [23 , with a toy model explaining the construction. But this 
is just a toy model. For a more realistic, but still not too unwieldy model 
of an algebraic implementation of the constraint, T> => p , let T>q be just a 
projector of rank one, | V E'}(\E , |. Use the spherical harmonic oscillator with its 
standard set of orbitals, ip n i m ar(r) = Yi m (r) exp(-j/r 2 /2) P n i(r) Xa Xr , where 
a = ± and r = p,n are spin and isospin labels, and, more important, the orbitals 
have the same exponential decay, modulated by polynomials P n i . Truncate this 
single particle basis at some maximum value of the number of hu; quanta. The 
truncation may differ for neutrons and protons, if necessary. Given a neutron 
and a proton numbers, N, Z, prepare an orthonormalized shell model basis of 
states $i for a configuration mixing, \& = ^ . (cj+i c'j ) $ j , to describe this N, Z 
nucleus, assumed to have a ground state spin zero. (Here the real and imaginary 
parts, Cj, c'j of the mixing coefficients have been explicited.) The $i's, naturally, 
are made of iV + Z orbitals (fulmar, with the necessary recoupling of angular 
momenta to induce total angular momentum 0. Given H, tabulate the matrix 
elements, i?y = making most often a sparse matrix. The energy is, 




0) 



Then the neutron density, 
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Pn(r) ^ ^ (4+n a f+n + 4-n a r-n) (Cj + I c'j), (10) 

ij 

is also a polynomial in terms of the real and imaginary parts c's, c"s. It is also a 
radial function, p{r). Moreover, it is a polynomial of r, multiplied by exp(— vr 2 ). 
Its order Af n is indeed finite, since the single particle basis was truncated at 
some harmonic oscillator major shell. Only N n of the Af n + 1 coefficients are 
independent, because the integral, J Q r 2 dr p n (r), equates the neutron number. 

Let CJ^ denote the independent coefficient multiplying r m in this "neutron 
density polynomial" . It can be stressed that the matrix element , \at a? n | $j) , 
(spin summation understood) is also the product of exp(— vr 2 ) by a polynomial 
of r. Hence, let such matrix elements be tabulated, and, in particular, let their 
coefficients C™ m of r m be tabulated. The "neutron density constraint" , in the 
Levy and Lieb sense pQ, is thus the set of algebraic relations, 

5> 4 - 1 4) ( Cj + 1 4) C? jm = , m = 1, M n . (11) 

The same considerations hold for the proton density constraint, leading to 
the set of polynomial constraints, 

^(c i -id i ){c i +it! j )C p ijm = CP m , m=l,...,Af p . (12) 

ij 

The last constraint to consider is the normalization of "P, 

E^ + (cD 2 = l, (13) 

i 

hence, the total number of constraints is, a priori, J\f c = M n + A/" p + 1. 

A priori also, the total number of real parameters Ci, c! i is twice the number 
of states $i, diminished by 1, because of the arbitrary phase of Accordingly, 
one of the c"s, for instance, can be frozen to zero. Actually except for clumsy 
choices of phases in the basis, or special Hamiltonians breaking time reversal 
symmetry, it is reasonable to assume that all c"s vanish. The model can be run 
with c's only, the number N of which is the dimension of the basis {'I'i}. 

In every realistic model one can think of, the number N of parameters is 
significantly larger than that, Af c , of constraints. Then, because of Eqs. (fTTj) . 
(|T2|) and (fT3)) , one can eliminate Af c among the coefficients Cj from Eq. ([9]). 
This leaves a "precursor" polynomial relation, 

pprUv, c?, c», cy, ci c^ p , c r v ) = o, (w) 

between the energy i], the coefficients CJ^ P of the density polynomials and the 
(N — M c ) "survivor" coefficients cf urv . Then the minimization of ij under the 
constraint of a given density pQ induces (N— M c ) equations, d7 , prec /dcf urv = 0. 
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These are used in turn to eliminate the survivors. Every step in the process is 
an easy manipulation of polynomials. One finally obtains a polynomial relation, 

TZ(r], C?, CJ, C n Mn , C?, Cf , CfrJ = 0, (15) 

which links the minimal energy to the parameters of p. Although not an open 
formula, this implicit equation, 1Z = 0, has all the needed virtues of a density 
functional for the calculation of r\ and its minimization in terms of the density. 

This detour via the density to avoid a direct diagonalization of the matrix 
{Hij} is of some interest if the matrix is sparse but too huge to be handled 
directly. It is not excluded, in particular, that those polynomials 1Z obtained 
from "small" models turn out to give reasonable approximations for "bigger" 
models. It is indeed likely that physics allows a natural cut-off in the order of 
the density polynomials, because high orders might mean oscillations that are 
too costly in energy to be acceptable for ground states. 

It is clear that this rigorous derivation raises several technical problems, 
such as, for instance, the best tactics to handle big polynomials in the con- 
text of many parameters, the comparison of distinct polynomials 1Z obtained 
from models with different choices for the basis, the search for extrapolations 
when calculations become too cumbersome. But there is no doubt that this 
constructive approach offers challenging open problems. 

An interesting side question arises: can this algebraic method also provide 
the KS potential? We have preliminary results which are optimistic, but they 
are only ... preliminary Another interesting question is, can this algebraic 
method be used for the cluster model? We have again preliminary results which 
are optimistic. But they are also ... too preliminary. 

Last, but not least, it can be noticed that this polynomial method can be 
used for other collective degrees of freedom than the density. Every quadrupole, 
octupole, etc, indeed boils down to a polynomial with respect to the configura- 
tion mixing coefficients. 

6 Discussion and Conclusion 

Five problems have been stated by this note. The most urgent seems to be 
the concavity question. Constrained minimizations in a landscape that is not 
concave can lead to absurd mistakes. 

Then comes the question of the radial theory The empirical success of the 
present three-dimensional approaches, EDF in particular, should not hide the 
fact that these contradict the rotational invariance of the Hamiltonian. The 
excuse that mean field theory can break symmetry is correct, but is not accept- 
able for a rigorous DFT. Density functional theory is not a mean field theory, 
despite the Kohn-Sham formulation. This KS formulation must be readapted 
to the scalar world demanded by the scalar nature of H. 

It might be that three-dimensional approaches can be justified within the 
third subject discussed by this note, a variational principle based upon variation 
after projection of good quantum numbers. But this is yet an open subject. 
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Unless we learn of further evidence, the fourth subject, about the square root 
of the density, does not seem to deserve priority. Still, it would relate various 
effective potentials, of which nuclear physics is an active consumer. 

The last subject, a constructive method via polynomials, stands out as spe- 
cial. It uses a highly non local parametrization of the density, that deviates from 
the (quasi) local tradition of the field. This untraditional parametrization of p, 
in the frame of a polynomial algebra, introduces a completely new zoology of 
polynomial DFs. Everything must be reinvented in this context, for the number 
of soluble models is huge, calculations can be heavy and it is not obvious how to 
take advantage of such rigorous solutions ... that might be sometimes obscure. 
But isn't this a stimulating situation? 
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